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ABSTRACT 

We report the detection, with a high level of confidence, of coherent outflows around voids found in the seventh data release of the 
Sloan Digital Sky Survey (SDSS DR7). In particular, we developed a robust < |cos 6\ > statistical test to quantify the strength of 
redshift-space distortions (RSD) associated with extended coherent velocity fields. We consistently find that the vector that joins void 
centers with galaxies that lie in shells around them is more likely to be perpendicular to the line-of-sight than parallel to it. This effect 
is clear evidence for the existence of outflows in the vicinity of voids. We show that the RSD exist for a wide range of void radius 
and shell thickness, but they are more evident in the largest voids in our sample. For instance, we find that the < |cos 8\ > for galaxies 
located in shells within 2/r'Mpc from the edge of voids larger than 15/?~'Mpc deviates 3.8 lo" from uniformity. The measurements 
presented in this work provide useful information to constrain cosmological parameters, in particular Q.„, and erg. 
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1. Introduction 

According to the cosmological principle, the clustering of galax- 
ies is statistically isotropic. In actual galaxy redshift surveys, 
anisotropies arise because distances to galaxies are obtained 
from redshifts (i.e. Hubble flow plus peculiar velocities), making 
the line-of-sight a preferred direction. These so-called redshift- 
space distortions (RSD) pose a problem to several large-scale 
structure (LSS) studies, especially those that rely on the 3-d 
positions of galaxies (see Hamilton 1998 for a comprehensive 
review of the effect). Yet, RSD provide valuable information 
about cosmological parameters, in particular the matter density 
parameter Q. m and the normalization of the amplitude o f den- 



sity fluctuations erg (e.g. Kaisedll987 ; Cole, Fisher, & Weinberg 
1994tlTinker. Weinberg. & Zhengl2006HPercival & White 2009; 



Ross et al] 1201 lb . The RSD studies have improved consider- 
ably since early galaxy surveys thanks to measurements made 
in recent large-scale surve ys both from the analysis of the 
galax y correlation functions (Peacock et al. 2001 ; Hawkin s et al 
2003 using the 2dFGRS and Zeh avi et alJl2005HOkumura et al 



2008; Cabre & Gaztanaga 2009 in the SDSS) and power spec- 
trum (Per cival et al.ll2004l in the 2dFGRS). These measurements 
have also been ex tended to higher redshift dGuzzo et al.L [2008; 
Blak e et al.L 120101) . However, it is still debated which the best 
RSD estimator is for extracting cosmological information (see 
Percival & White 2009 for a review of this topic). 

A particularly interesting case for observing and quantify- 
ing RSD associated with extended coherent velocity fields may 
be provided by large cosmic voids in the distribution of galax- 
ies. These regions in the LSS have, on average, a very low en- 
closed matter density (5 = p(< R void ) - p mecm )lp m ean - -0.9, 



e.g. lGottlober et al.ll2003r.lPatiri et al.ll2006al) . which generates a 
well-defined spherically averaged radial outflow in the matter 
(and galaxy) velocity field (e.g. [P adilla, Ceccarelli, & Lambas 
l2005UCeccarelli et alJ l2006t iPatiri et alJ l2006bl) . As a conse- 
quence, the average density contours around large voids that 
are nearly spherical in real space will appear as ellipsoids elon- 
gated along the line-of-sight in redshift space. That is, galaxies 
mapped in actual redshift surveys that are located in spherical 
shells around voids are more likely to have a position vector per- 
pendicular to the line-of-sight than parallel to it. Quantifying this 
effect provides a powerful statistic to assess the strength of RSD. 
Moreover, a key advantage of studying RSD around large voids 
is that the strength of the effect can be computed efficiently for 
a given power spectrum using an existing theor etical framework 
(IPatiri et al.Ll2006aHBeTancort-Riio et al.Ll2009h . which is essen- 
tial to establish strong constraints to cosmological parameters. 

In this paper, we use a sample of galaxies located in the 
vicinity of large voids drawn from the SDSS DR7 to study 
redshift-space distortions. To quantify RSD, we employed a 
modified form of a particularly simple, yet robust method: the 
< |cos 9\ > statistics, i.e. the averaged absolute value of the 
cosine of the angle, 8, between the position vector of galaxies 
within shells around voids and the line-of-sight to the center of 
their corresponding void. This provides a clear, powerful statistic 
to assess the strength of RSD. 

The paper is organized as follows. In section 2 we present the 
data and sample selection. In section 3 we describe the statistical 
test and show our results, and finally in section 4 we discuss the 
results and future work. 
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Table 1. Complete statistics of large voids found in our volume- 
limited galaxy sample extracted from the SDSS DR7. Afy de- 
notes the number of voids larger than the given radius. 



Radius 


N v 


(/r'Mpc) 




10 


1162 


11 


829 


12 


464 


13 


284 


14 


141 


15 


76 


16 


34 


17 


20 



2. Data and sample selection 

We used the NYU V alue-Added Galaxy Catalog (NYU-VAGCfl 
iBlanton et al.ll2005l) associated to the Sloan Digital Sky Survey 
Data Release Seven (SDSS DR7, Abazaiia n et al. Il2009h . which 
is the final data release for the SDSS-II. The full catalog includes 
redshifts of ~ 700, 000 main sample galaxies covering ~ 8000 
deg 2 on the sky. From this catalog, we extracted a volume- 
limited sample defined by a r-band absolute magnitude thresh- 
old M r - 5 log h - -20.3 and a maximum redshift of z = 0.13, 
which resulted in a sample with a total of 162,076 galaxies. Note 
that all galaxies are fc-corrected t o redshift z = 0-1 (a pprox. the 
median redshift of the survey, see lBlanton et al.ll2003l) . 

We assumed a flat ACDM cosmology with parameters con- 
sistent with the latest WMAP-7 results djarosik e t al. 201 1). We 
converted galaxy reds hifts to prop er comoving distances using 
the standard equation dHoggll 19991) 



H Jo V^m(l+Z') 3 +0\ 



(1) 



with Q m = 0.267 and Qa = 0.733. We write the Hubble constant 
Hq in terms of the dimensionless Hubble parameter h, i.e. Hq - 
100/l We checked that galaxies are uniformly distributed within 
the sample applying a standard V max test. 

In this work, we define voids as maximal spheres empty of 
galaxies brighter than the limiting magnitude of our galaxy sam- 
ple defined above. To search for voids in the gala xy sample, 
we used the HB void finder. We refer the reader to iPatiri et aD 
(2006b) for details about the void-finding methodology and al- 
gorithm. In Table [TJ we present the statistics of voids found in 
our volume-limited galaxy sample. 

Here, we are interested in the galaxies located in the bound- 
aries of voids. In particular, we searched for all galaxies located 
in shells of a given width around the edges of the voids. Then, 
we computed the angle, 9, between the position vector of the 
galaxy (with respect to the center of the void) and the position 
of the center of the void as follows 



s • r 

WW 



(2) 



where s is the position vector of the galaxy with respect to the 
center of the void, and r is the position vector of the center of 
the void with respect to the observer. 



3. Statistical test and results 

In principle, variations of standard methods, such 
as the void -galaxy correlation function (see e.g. 
Padilla, Ceccarelli, & Lambas 20051) can be used to detect 
and quantify RSD around voids. However, in the largest cur- 
rently available galaxy survey such as the SDSS, the number of 
voids and galaxies in their vicinity are not sufficient to obtain 
sharp density contours and consequently reliable statistics 
necessary to extract inferences for the relevant cosmological 
parameters. As mentioned above, we developed a modified 
version of the widely used < |cos 9\ > statistical test. 

The standard way to construct a < |cos 9\ > test would be 
to simply calculate the average of |cos 9\ over all galaxies in the 
sample located in the relevant shells of a given width around 
the edges of the voids. However, because of the correlations 
between the galaxies within a shell of a given void, the corre- 
sponding cosine values are not independent random variables. 
This hinders assessing the statistical significance of the test. For 
instance, those correlations may cause the rms to increase con- 
siderably from the standard 1 / ^j\2N g rms (where N g is the total 
number of galaxies in the void shells) obtained if galaxies were 
uncorrelated. 

We constructed a modified version of the standard average 
of the |cos 9\ test in the following way 



1 Ny 

<<; |cos 9\ »= — V < |cos 9\ >, 



(3) 



where < |cos 9\ >,■ is the average of |cos 9\ over all galaxies in 
the relevant shell of the j-th void and Ny is the total number of 
voids in our sample. 

In this procedure, galaxies located in the shell of a given void 
are distributed independently for different voids, so they can be 
treated effectively as an independent random variable. To esti- 
mate the rms of our statistic we used the independence of all 
terms of the sum in the r.h.s. of equation [3] so the variance of 
our statistic is given by 



rms(«: |cos 9\ ») = — x 
Ny 



(2S<|cos0|> 2 ) 



N v 



— <sC \cos9 » 



(4) 
(5) 



Once we obtained the rms we computed the statistical sig- 
nificance (i.e. giving the confidence level at which the null- 
hypothesis may be rejected) written as the number of sigmas, 
No-, in the following way 



N^ = 



(« |cos 6»| » -0.5) 
rms(<z: |cos 9\ ») 



(6) 



http://sdss.physics.nyu.edu/vagc 



It is worth noting that this is a conservative approach. In prin- 
ciple, it should be possible to construct a more efficient statistical 
test. This task, however, is not simple to carry out since the sam- 
pling errors of that estimate are difficult to treat. We will address 
these issues in detail in a future work. 

In Table 2 we present the results for the <K |cos 9\ » statistic 
and its respective significance. We show the results obtained for 
several shell sizes and void radii. In absence of redshift distor- 
tions, i.e. when the redshift space coincides with real space, the 
average value of the absolute value of cos 9 must be 0.5. 
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Table 2. Summary of our results. 









Radius 
(h-'Mpc) 


<k |cos 0\ » 


rms 


N LT 
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N s 


< |cos e\ > r 


rms r 


Shell width= 


=2.o/r 


'Mpc 
























10 

12 
15 


0.4886 
0.4934 
0.4509 


3.758E-03 
6.059E-03 
1.288E-02 


3.02 
1.09 
3.81 


1.047+U.UJ6 

1.027 + ^ 

1 T1 0+0:052 


16794 
6926 
1120 


0.4991 
0.4973 
0.4940 


4.7810E-03 
6.8347E-03 
1.6674E-02 


Shell width= 


-4M' 


'Mpc 
























10 

12 
15 


0.4919 
0.4935 
0.4613 


3.001E-03 
4.856E-03 
1.123E-02 


2.70 
1.34 
3.44 


1 033 +(u,li 
1 026 +!llM 
1.168_ nf)S4 


41837 
17297 
2791 


0.4990 
0.4964 
0.4903 


3.4634E-03 
5.7295E-03 
1.3100E-02 


Shell width= 


=6.0/?- 


'Mpc 
























10 

12 
15 


0.4949 
0.4933 
0.4651 


2.593E-03 
4.182E-03 
1.005E-02 


1.96 
1.60 

3.47 


1.021 +ou " 
1 027+^'^ 

'• ui ' -0.017 

1 150 +0 - 045 

i.ijo 0047 


76887 
32062 
5243 


0.4999 
0.4986 
0.4899 


2.8861E-03 
4.5821E-03 
9.5722E-03 


Shell width= 


=8.0r 


'Mpc 
























10 

12 
15 


0.4966 
0.4927 
0.4711 


2.331E-03 
3.749E-03 
8.666E-03 


1.46 
1.96 

3.33 


u -MUD 

1 030 +00 ' 5 
1 191+0.038 


121031 
51286 
8520 


0.5000 
0.4977 
0.4900 


2.5342E-03 
4.2074E-03 
9.9732E-03 



Notes. <K |cos 8\ » is the average over all voids of the average of |cos G\ of the galaxies in the shells of the void sample larger than the given radius. 
N a is the number of sigmas for the statistics. N g is the total number of galaxies in the given shell, p is the anisotropy parameter. < |cos 6\ >,- shows 
the value of |cos G\ for galaxies in the shells of randomly placed spheres and rms,. is its root mean square calculated over 100 realizations. See text 
for details. 



We double-checked this by placing random spheres over the 
galaxy sample and calculating the angles for the galaxies found 
in the shells of the random spheres following equation [3] In 
particular, we made 100 different realizations of Ny random 
spheres, where Nv is the number of actual voids larger than a 
given radius. We then calculated the standard average and rms 
of the 100 realizations. We show these results also in Table 2 
(columns 6 and 7). Note the similarity between the rms values 
shown in columns 3 and 7. This reflects the fact that the corre- 
sponding settings, i.e. the number and correlations of galaxies 
in the shells of voids and the random spheres, are rather similar. 
Also, because the values of < |cos 9\ > for the random spheres 
are 0.5 (within the rms) possible biases caused by border effects 
or by not using the "correct" cosmological parameters in equa- 
tionQ]can be excluded. 

For the wide range of void radii and shell sizes considered, 
the statistic generally shows values of «: |cos 9\ » consistently 
lower than 0.5. This departure of the distribution of «: |cos 9\ » 
from uniformity is a clear signal of RSD generated by the under- 
density associated with the voids. The effect is a unique func- 
tion of the value of the enclosed fractional density contrast (6) 
within the average radius of the galaxies within the shell, and 
does not depend directly on the value of that radius. Moreover, 
the largest deviations appear for the largest voids in our sample. 
This is not surprising because the average underde nsity within 
voids increases with the radius of the void (see e.g. IPatiri et alJ 
2006a), and as mentioned above, if the enclosed density within 
voids decreases, the radial outflows will be stronger. 

On the other hand, the effect is more pronounced for the 
smaller shell sizes. This trend can be interpreted by checking 
the profiles of the en closed density inside and around the voids . 
Several works (e.g. [Gottlober et al.l 120031 IPatiri et al.l r2006bt 



ICeccarelli et alj|2006l) showed that the void density profiles are 
monotonically increasing. Consequently, for a given sample of 
void radii the effect must be monotonically decreasing (safe for 
statistical fluctuations) as the width of the shell increases. We 
obtain, in general, considerable values for No- both for large and 



smaller voids, but the largest sigma is found for galaxies in shells 
of 2.0/r'Mpc of voids larger than 15.0/r'Mpc. 

To show our results in a more intuitive way, we use an 
anisotropy parameter, p, defined by 



1 



<K |cos 6*| » 



1. 



(7) 



Within this parameterization, p values lower than 1 mean that 
the vectors s and r tend to be parallel, while for p values higher 
than 1 the vectors tend to be perpendicular, p = 1 means com- 
patible with isotropy. We also show the values of p for the cases 
considered above in Table 2, showing more patently the trends 
seen with the <K |cos 9\ » test. 

It is also interesting to analyze the probability distribution of 
the |cos 9\ values, defined as P(|cos 9\) (hereafter P, for brevity). 
If there are no distortions, P must be flat, equal to 1 . In figure 
Q] we show as an example the P distribution for voids larger 
than lS/r'Mpc and galaxies within shells of 2/?T 1 Mpc (filled 
circles). We can clearly see that the distribution has a declin- 
ing trend, incompatible with isotropy. In rece nt theoretical works 
dCuesta et al.ll2008l:lBetancort-Riio & Truiilloll2009h . it has been 
shown that a general probability distribution for quadrupo lar dis- 
tortions is adequately described by the following equation: 



P(|cos 9\) d\cos 9\ = 



p d\cos 9\ 



(l + (/? 2 -l)|cos<9| 2 ) 3 / 2 ' 



(8) 



where p is the anisotropy parameter defined in equation (|7). 
We also show in figure Q] the predicted probability distribution 
(dashed line) for the case studied above, i.e. with p - 1.218, 
which agrees well with the observational results. 

4. Discussion and conclusions 

Our results show clear evidence for the existence of outflows 
in the vicinity of large cosmic voids found in the SDSS DR7. 
Because «: |cos 9\ » is smaller than 0.5 for the wide range of 
void radii and shell sizes computed in our sample of galaxies 
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Fig. 1. Probability distribution of |cos 9\ for voids larger than 
15/i^'Mpc and galaxies located in shells of 2/z~ 1 Mpc (filled cir- 
cles). The full flat line indicates the distribution for isotropy. The 
dashed line is a theoretical prediction, see text for details. 



trum, which makes large density fluctuations more likely, which 
in turn increases the value of |<5|. 
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and voids, the vector that joins the galaxies located in the shells 
around voids with the centers of those voids is more likely to be 
perpendicular to the line-of-sight than in a random distribution. 
This effect is understood by means of redshift-space distortions: 
in redshift space, our spherical shells intersect higher density iso- 
contours in the direction perpendicular to the line-of-sight (recall 
that a sphere in real space is elongated along the line-of-sight due 
to redshift distortions). Furthermore, the effect increases for void 
with larger radius (for a given shell width), while for a fixed void 
radius it decreases for a broader shell width. 

We may interpret this in terms of the relationship between 
the deformation of the density contour, T, and the underlying 
density contrast, 6, which can be written as follows 



T cc f(6)H 



dlnD 
d lna ' 



(9) 



where D and a are the growth and scale factors, respectively, 
Hq is the Hubble constant and f(5) is a certain function of 
the enclosed fractional den s ity fluctuation 5 within the radius 
R (ISheth & Tormenl 120021: iBetancort-Rijo et all 120061) . This 
function monotonically decreases with the absolute value of 6. 
Therefore, for a fixed shell width, the effect is more significant 
for larger voids (since they have higher values of 5). For a fixed 
void radius, the effect is roughly proportional to f(6), with 6 
equal to the enclosed fractional density fluctuation within the 
mean radius of the shell. Thus, as the shell width increases, the 
effect decreases since the average value of \6\ will be lower. 

In a forthcoming paper, we will explore in detail the model 
predictions for the magnitude of the effect of redshift distortions 
around voids as a function of cosmological paramete rs. For this 
we will use the analyti cal framework presented in iPatiri et all 
(l2006ah and extended in lBetancort-Riio et all (120091) . In particu- 
lar, we will study the dependence of the p parameter with Q,„ 
and cr g . The dependence with Q m enters through d InD/d Ina 
(roughly proportional to QjJ; 6 ) and also via F = Q, m h, which 
characterizes the shape of the power spectrum. The dependence 
on erg is the result of the higher amplitude of the power spec- 
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